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PREFACE 

The contents of this report formed the basis for a lecture given 
by one of the authors at the Freshwater Institute in March, 1977. The 
work described sunmarizes some important aspects of the research carried 
out by the Limnology and Plankton Taxonomy Units since their inception; 
thus it was decided that it v/ould be useful to release ttn's information 
in a report form. 
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ABST".ACT 



Freshwater managenent lias suffcrec from a lack of sinple predictive 
models. It is desirable that such F.ociels be easily manipulated Ly the 
lake manager and offer hin credibility end confidence v;fien it comes to 
making decisions. Improved modelling methodology for phosphorus is 
presented along vnth further evidence to support the hypothesis that 
phosphorus is the key element in tiie biogeochemistry of north temperate 
lakes. It is stressed that the lake and its watershed are the unit under 
consideration wlien employing such n-odels. Once phospliorus concentration 
is lietermint'd, other water quality parameters such as algal biomass and 
nypolimnetic oxyoen depletion rate can be predicted. Ultimately, one 
wnuld like to employ a master variable such as phosphorus concentration 
from wnich to extrapolate many ottier aquatic ecosystem parameters but 
tlie complexity of dynamic natural systers limits this goal. Further 
discussion emphasizes the complexity of nutrient-phytoplankton relationships 
in lakes and that, as a result, lake mananenent techniques shoulc' still 
be practiced witfi caution. 



DISCUSSION ANP SLTIMARY 

One of the important questions that our group in the Ministry of 
the Environment is trying to answer is "How much developrtont, in terms 
of houses, cottages, condominiums or campgrounds etc., can a piven body 
of water withstand before degradation of water quality occurs beyond an 
acceptable level or, conversely, how much will any specified human activity 
quantitatively effect the water quality of any body of water associated 
with it?", and to describe our present attempts to answer this question, 
we have to review a number of points. 

It is now widely recognized tliat P is the most important element 
governing the trophic status of most lakes. We have results froni three 
studies which reinforce this hypothesis very clearly (Nicholls et al. 
1977, Dillon et al. 1978, Scheider and Dillon 1976). 

The first example is Lake Erie's Western Basin, The accelerated 
eutrophication of Lake Erie over the period 1919 to 1965 has been well 
documented by Davis (1964), Verduin (1964) and Hohn (1969) who summarized 
long-term phytoplankton records for the lake. By the late 1950's and 
early 1970' s, attempts were already underway to halt the eutrophication of 
Lake Erie by decreasing the phosphorus load to the lake. It had been 
estimated that phosphorus of detergent origin contributed about 40% of 
ttie total phosphorus loading to Lake Erie, prior to legislative action by 
both Canada and the U.S.A. between 1971 and 1973 to limit thk.> pliosphorus 
content of detergents to 8.7% (as P, by weight). At weekly intervals 
since 1967, phytoplankton analyses riave been performed on water collected 
at the Union Water Treatment Plant in the Western Basin of Lake Erie. Since 
1971, all classes of algae have declined sharply and have accompanied 
decreased concentrations of total phosphorus in the Western Basin near-shore 
area (samples collected also by the Union Hater Treatment Plant) and 
decreased phosphorus loading from the Detroit River (Hiciiigan Jepartment 
of Natural Resources data). 

Since the 1967-1970 period, when average phytoplankton density was 
5010 A.S.U. ml"' phytoplankton density had fallen about 42?: to 2900 A.S.U. 
nil~^ by 1975. There was no further change in 1976 and by 1977 algal density 
had increased slightly (Fig. 1). Decreased densities of diatoms account 
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Fig. 1. Average annual phytoplankton density {Area! Standard Units/inl) 
and total phosphorus concentration in near-shore western Lake 
Erie, determined from samples collected approximately weekly 
at the Union Water Treatment Plant at Kingsville, Ontario 
{updated from NichoUs et al^. 1977). 



for most of the decrease in total phytoplankton, since this group comprised 
65-80% of the total during the long-term study. These findings should 
offer encouragement for the "reclamation" of other highly eutrophic inshore 
waters of the Great Lakes such as Lake Ontario's Bay of Quinte, Lake 
Huron's Saginaw Bay and Lake Michigan's Green Bay, given similar controls 
on phosphorus loading. 

The second example is Gravenhurst Bay (Fig. 2) which lies on the 
Precambrian Shield and which has been extensively studied from 1969-77 
(Michalski and Nicholls 1975; Michalski et al. 1975, Dillon et aji^. 1978). 
In late 1971, phosphorus removal was initiated at the town's sewage 
treatment plant as a test case for nutrient diversion effects upon a lake. 
From Table 1 it can be seen that by 1975 removal was effectual in 
reducing P concentration in the Bay to 45* of the pre-removal average, 
chlorophyll a_ concentration to 51 "J, the quantity of phytoplankton to 48« 
and the areal hypolimnetic oxygen demand to 50%. From Fig. 3 it is clear 
that the areal hypolimnetic oxygen demand in Gravenhurst Bay was much 
greater in 1969 than 1975. 

By 1975, the Bay had shown good overall improvement, however, in 
June of 1976, a massive bloom of Aphanizomenon occurred. Apparently this 
was the result of the failure of P removal at a second sewage treatment 
plant which was brought into operation in October 1975 tc accommodate 
the expansion of Gravenhurst. About 50% of the town's wastes were directed 
through the new sewage treatment plant which was only operating at 20-30% 
P removal efficiency instead of the expected 80-90%. The failure continued 
from October, 1975 to April, 1976 at which time it was rectified. The results 
were severe (Table 2): chlorophyll a_ levels and areal hypolimnetic oxygen 
demand were as high and phytoplankton quantities higher than during the 
pre-removal period. P concentrations and Secchi disc values were only 
slightly worse than in 1975 but these mean values do not adequately represent 
the radical changes within the Bay. For instance, even though Secchi disc 
was the lowest ever recorded in the early summer of 1976 (Fig, 4), by fall 
it was the best. The same pattern showed up in chlorophyll a_ levels (Fig, 5) 
and P concentrations (Fig. 6). Similarly, the phytoplankton community 




10 km. 



79: 
24 



Fig. 2 Location and morphometric map of Gravenhurst Bay. Sewage 
treatment plants are marked with a triangle, sampling 
station with a star. 



Table 1: Mean total phosphorus, chlorophyll a, Secchi disc depth (S.D.), 
quantity of phytoplankton and areal hypolimnetic oxygen deficit 
(AHOD) in Gravenhurst Bay for three years prior and four years 
following P removal (from Dillon^ al_. 1978). 

Year [P] [chl a] S.D. Phyto. AHOD 

(mg.m"M (mg-m"^) (m) (ASU-ml'M (mg.cm~^d-M 

1969 42 10.6 2.6 1780 0.053 

1970 39 5.1 3.1 1620 0.047 

1971 52 13.8 1.9 3620 



1972 35 8.1 3.1 

1973 33 6.9 3.2 

1974 25 5.0 2.7 

1975 20 5.0 3.9 



Mean 44 9.8 2.5 2340 0.050 
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Fig. 3. Mass of oxygen in hypoliinnion of Gravenhurst Bay as a function of time before (1969) and after (1975) P 
removal. Slope of the line is the areal hypolimnetic oxygen deficit rate. 
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Table 2: Comparison of trophic status parameters (defined in Table 1) 
measured in Gravenhurst Bay before P removal with those 
measured in 1975 and 1976 (in part, from Dillon et al. 197s). 

Year [P] (chl a] S.D. Phyto. AHOD 

(mg.tn-^) (mg.m"^) (m) (ASU.ml'M (mg.cm-^.d"M 

1969-71 44 9.8 2.5 2340 0.050 



1975 20 5.0 3.9 1120 0.025 



1976 27 10.6 3.7 4980 0.046 
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Fig. 4 Secchi disc depths in Gravenhurst Bay, 1975 and 1976. 
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Fig. 5 Chlorophyll a^ concentration in the euphotic zone of Gravenhurst 
Bay (measured as twice the Secchi depth) from 1974-76. 
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Fig. 6 Total phosphorus concentration in the euphotic zone of 
Gravenhurst Bay, 1974-76. 
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changed from one dominated by Aphanizomenon in June and July to one 
resembling oligotrophic Precambrian Shield lakes with low numbers of 
diatoms, motile Chrysophyceae and Cryptophyceae (Fig. 7). 

Although changes since 1971 in phytoplankton biomass in Gravenhurst 
Bay have been dramatic and are relevant to many aspects of water quality 
(e.g. dissolved oxygen depletion - see Fig. 8), changes in phytoplankton 
composition are equally significant. 

Declines in the Cryptophyceae and Chlorophyceae and increases in 
the Chrysophyceae have occurred independently of the large fluctuations 
of the Cyanophyceae or blue-green algae (Fig. 9). The blue-green algae, 
especially the bloom-forming genera, Aphanizomenon , Anabaena and Microcystis 
seem best able to respond quickly to changes in nutrient loading and were 
found in high densities during the summers of 1971 and 1976 as a consequence 
of high nutrient concentrations. This independence of the blue-green algae, 
especially the bloom-forming genera, from the rest of the algal flora is 
best illustrated by Fig. 9, which shows that the group is an important 
component of the plankton only when present at high density. During years 
when blue-green density is low, other algae are much more prominent in the 
plankton. 

The blue-green algae have responded most to changes in total 
phosphorus loading and concentration; however, it would seem important to 
separate the bloom- formers ( Aphanizomenon , Anabaena , Microcystis ) from the 
"non-bloomers" ( Aphanothece , Aphanocapsa , Merismopedia and Chroococcus ) 
since the response of these two groups to the phosphorus control programme 
has been quite different. By 1975, the average blue-green "bloomer" 
biomass had declined approximately 78% to 170 A.S.U.ml'^ from a pre- 
phosphorus control period average of 770 A.S.U.ml".' In contrast, the 
"non-bloomers" had increased from an average of 40 A.S.U. ml"^ during the 
1969-1971 period to 270 A.S.U. ml"' during 1975 (Fig. 10). 

In other lakes in Ontario, certain Chrysophyceae (mainly species 
of Dinobryon , Uroglena , Mallomonas and Synura) are found in the metal imnion 
or upper hypol imnion of lakes along with several cyanophycean "non-bloomers" 
(e.g. Fig. 11 and unpub. data). In view of the similar niche occupied 
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Total density (Areal Standard Units/ml and composition of phytoplankton in Gravenhurst Bay during 1976 
As a result of sewage treatment plant malfunction, the phytoplankton during the June-July period are 
indicative of extreme eutrophy. The rapid decline in algal density and alteration in compositio 
correction of sewage treatment plant malfunction, created typical oligotrophic characteristics ^" 
plankton during the late suinner and fall period. 
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ALGAL BIOMASS (XIO^SU) 



Fig. 8. Correlation between areal hypolimnetic oxygen deficit (AHOD) 

and average euphotic zone algal density (Areal Standard Units/ml) 
during a five year period in Gravenhurst Bay of Lake Muskoka. 
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Fig. 9. Changes in average density of several important classes of algae 
in Gravenhurst Bay from 1959 through 1976, 
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Fig. 10. Relationships between average total blue-green algal density ( Cyanophyceae ) and blue-green "bloomer:" (i.e. 
mainly Anabaena , Aphanizomenon and Microcystis spp. ) as a percentage of the total phytoplankton. 



5- 



£ 10 
O 



15- 



Mallomonas caudata 







60 



100 



(a) 



(b) 




T 




50 



(0 



Oinobryon bavaricum 







200 



"T — 
400 



Cell Volume (mmVm^) 



Fig. n 



Vertical distribution of some common members of the Chrysophyceae in (a) Buckskin Lake (78° 12' W 45° 57' wl 
(b) Tadenac Lake (79° 57' W, 45° 03' N) and (c) Buchanan Lake (78° 50' W. 45° 17' N). ^' 
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by both groups, it is not surprising that their responses to the cluiiuie 
in nutrient input to Gravenhurst Bay would be similar. By 197^, the 
Chrysophyceae in the plankton of Gravenhurst Bay had increased to an ice- 
free period average of 120 A.S.U.ml"^ from a pre-phosphorus control period 
average of only 16 A.S.U.ml"' (Fig. 12). Both the Chrysophyceae and the 
blue-green "non-bloomers" would appear to be particularly well adapted 
to low light and temperatures. Increases in biomass of these fonns could 
only occur after significant declines in total phytoplankton (mainly 
diatoms and bloom- forming blue-green algae) in Gravenhurst Bay had allowed 
increased light penetration to the deeper waters, thereby allowing growth 
of forms more suited to these depths. 



In years prior to 1972, there was little inorganic nitrogen in 
Gravenhurst Bay during the summer period (Fig.13). Inorganic nitrogen built 
up as production declined while total nitrogen remained relatively constant 
over the entire season. The inorganic nitrogen disappeared in the following 
vernal period as it was incorporated into organic matter. In 1975, a year 
of P removal, the depletion of inorganic nitrogen did not occur over the 
summer period, indicating that it was not used by phytoplankton. This 
pattern was, however, reversed in 1976 after the sewage treatment plant 
failure and inorganic nitrogen distribution resembled pre-removal years 
(see also Table 3). Thus, response time to an increased input of P was 
very quick. Complete response to reduced input was also quick and appears 
to have taken about two months as can be seen from Fig. 14. The area! 
hypolimnetic oxygen demand plot can be broken down into two periods for 
1976; up to the end of June (0.088 mg O2 cm"^ d~^) and fron July to August 
(0.030 mg cm~^ d"'). This demonstrates the rapid response time generated 
by altering P inputs to lakes. Gilbertson, Dobson and Lee (1972) have 
shown a similar quick response time to nutrient loading and areal hypolimnetic 
oxygen demand in Lake Erie. 

In summary, improvement in Gravenhurst Bay was proportional to the 
P reduction. Thus P removal can be critical in treatment of effluents to 
some waterways. 

Middle Lake (Fig. 15) is a third case which tias inportant qualitative 
as well as quantitative results. We have studied some of the acidic lakes 
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Fig. 12. Changes of some constituents of the Gravenhurst Bay phytoplankton 
following controls on phosphorus loading from the local sewage 
treatment plant. 
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Table 3: Changes in total phosphorus content and the partitioning of 
nitrogen between inorganic and organic forms in the euphotic 
zone of Gravenhurst Bay , 1969-76. 



Year 


[P] 


[TIN] 


[TON] 


TN/TP 




(mg.m"^) 


(mg.m"^) 


(mg-m"^) 




1969 


42 


65 


384 


11.1 


1970 


39 


77 


427 


13.9 


1971 


52 


88 


527 


14.7 


1972 


35 


144 


383 


17.8 


1973 


33 


163 


399 


18.3 


1974 


25 


188 


356 


24.8 


1975 


20 


257 


320 


35.2 



1976 
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29 
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Fig. 13 Total inorganic nitrogen (TIN = NH3 + NO2 + NO3 ; solid line) 
and Total nitrogen (TN = TIN + organic nitrogen ; dotted line) 
in Gravenhurst Bay prior to phosphorus removal (1970), 1 and 4 
years after removal was initiated (1972 and 1975 respectively), 
and following a short-temi increase in phosphorus load (1976). 
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Fig. 14. Mass of oxygen in hypoiimnion of Gravenhurst Bay in 1976 as a function of time. Rate of decrease was high 
in May and June, much slower in July and August. 
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Fig. 15 Location of Middle Lake in the Sudbury region. 
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near Sudbury now for the last four years with the intetit of understandino 
the biological and chemical responses of these lakes to various chenical 
treatments, including decreased and increased pH as v/el 1 as nutrient 
additions. In 1973 fliddle Lake had a pfi of 4.2 which was raised to 
approximately 7 by addition of Ca(0fl)2 and CaCOj (Fig. 16). Middle Lake 
was fertilized in 1975 with P alone and this raised the P concentration 
from something near the detection limit of 2 yg 1~^ to 9 yg 1"^. This 
new 'high' post-fertilized level is still very lov/ when coripared to P con- 
centrations in oligotrophic Precambrian lakes. 

Although there was a quantitative biological response due to the 
fertilization of Middle Lake (e.g. ciilorophyll a_ and phytoplankton biomass 
increased) (Table 4), the significant change was qualitative. The phyto- 
plankton population in 1973 at pH 4.2 was dominated by cryptophytes , 
dinoflagellates, chrysophytes and chlorophytes (Fig. 17) which is typical 
of acid lakes in Sweden (Hornstrom et al_. 1973) and Ontario (Yan et al . 
1977). In 1974, at a pti of 6.5 to 7.0, the algal composition was altered 
drastically so that chrysophytes now dominated the phytoplankton community 
with diatoms appearing at different time intervals (Fig. 17), a seasonal 
pattern typical of circumneutral oligotrophic Precambrian lakes. Following 
P additions in late June of 1975, the community again drastically changed 
with blue-green algae completely dominating. Also during this time pH 
never exceeded 7, which is contrary to the theory of Shapiro 1973 that 
low pH excludes blue-green algae in lakes. Thus, this experiment clearly 
indicated the importance of P, even at what are normally considered to be 
•^ery low concentrations. 

Preliminary data from another acid lake, Mountaintop Lake, also 
in the Sudbury area (pH = 4.2), suggest that unless reserves of inorganic 
carbon are increased artificially as they were in fliddle Lake by neutra- 
lization with Ca(0H)2 and CaCOj , the expected yield of phytoplankton from 
additions of phosphorus are not achieved as quickly (nitrogen as well as 
phosphorus has been added to Mountaintop Lake). An average concentration 
of phosphorus (supplied as HjPO^) of about 60 yg P 1"^ was achieved by 
fertilization. However, after the first year, only about 15% of the 
expected yield of 9-10 mm^ l'^ phytoplankton was produced. IJe suspect that 
this lake is carbon limited under a regime of such higti rates of nitrogen 
and phosphorus additions. 
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Fig. 16. pH of Middle and Hannah Lakes, 1973-75. Triangles denote application of Ca(0H)2 and CaCO; 
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Table 4: Demonstration of the effects of fertilization on Middle Lake 
(from Scheider and Dillon 1976). 



Pre- Pre- 

fertilization fertilization 



Total phosphorus ^ r « -, 

(mg.m-3) ^-^ ^•" 



chlorophyll a , , 

(mg.m-3) " ^'^ 



Phytoplankton -..^ 

(A.S.U.-ml-M '^-^ 



Benthos 

{ n umbers -m"^) 



2v4 



703 



Zooplankton q g^ 2.32 

(numbers- 1 } 



96 138 



MIDDLE LAKE 




Fig. 17. Phytoplankton biomass (mg-l"^) and composition by major 
classes in Middle Lake, 1973-76. 
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It IS worth emphasizing that Ontario Lakes shov <\ s.-tv'\^ defined fe>la- 
tionship betv/een the average biomass of phytoplankton during the ice-free 
period and the average total phosphorus concentration (e.g. see Fig. 45). 
Given little year-to-year change in nutrient concentration, seasonal 
distribution of phytoplankton would also seem to be remarkably constant 
from year-to-year {Fig. 18). For many lake locations v/here we have several 
years of data, it is apparent that little year-to-year change in phosphorus 
concentration results in little change in average phytoplankton biomass 
(Fig. 19). In contrast, in those lakes which do experience a change in 
phosphorus concentration resulting either from lake fertilization or from 
improved waste treatment and decreased phosptiorus loading, the resultant 
change in phytoplankton yield can be predicted from the empirical data 
(Fig. 19). 

The next point to be made is that the parameter of in'portance is 

the P concentration in the lake, not the P loading. Vollenweider's 

original plot (1969) and its subsequent alteration (Vollenweider 1975, 1976) 

to L versus z/tw (Fig. 20) both have lines separating lakes into 

oligotrophic, mesotrophic and eutrophic categories based on a subjective 

assessment of lake conditions. The second plot, or one of its later 

variations, has definitely yielded a more realistic representation by 

including water replenishment time (x,,) or its inverse, flushing rate (p) 

w 

as shown in a study of Cameron Lake and Four Mile Lake in southern Ontario 
(Table 5) (Dillon 1975). Both lakes were of similar morphometry but 
Cameron Lake had an areal P loading (L) twenty times that of Four Mile Lake. 
The P and chlorophyll a^ concentrations and the Secchi disc depths, however, 
were virtually identical. The reason for the difference was that a high 
flushing rate in Cameron Lake offset the high loading rate. Thus, not only 
does the L versus z/t plot rationally give a better fit to empirical data, 

W 

as illustrated by the above case, but the lines separating the lake classes 
have slopes with appropriate units of mg-m'^-yr'^m-yr"' or mg-m"^. 
Therefore, both the rate and areal concepts are removed from consideration 
and lakes are separated according to some simple measure of their P 
concentration. 
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Fig. 18. Seasonal changes in total phytoplankton density at the outflow of Lake Ontario (Brockville) over a 
10 year period. 
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Fig. 19. Relationship between average total phosphorus concentration and 
average phytoplankton biomass (mm /I) at several lake locations 
in Ontario. Keyed as follows: E, western Lake Erie; 
Gr, Gravenhurst Bay of Lake Muskoka; M, Middle Lake near Sudbury 
(fertilized with H^P04); a, 5 years of data from the Bay of 
Quinte off Picton; b, 5 years of data from the Bay of Quinte 
off Conway; c, 5 years of data from Clearwater Lake near Sudbury. 
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Phosphorus load vs hydraulic load (qj). Lines designating "permissible" and "dangerous" loading levels as 
a function of qs are denoted. For further discussions, see Vollenweider (1975, 1976). 
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Table 5: Demonstration of the modifying effect of flushing rate on 

phosphorus loading. Two lakes of similar morphometry, one with 
a phosphorus load twenty times greater than the other, have 
similar lake phosphorus concentrations because of differences 
in flushing rates. 

Parameter Cameron Four Mile 

Aq (km^) 12.7 7.7 



(m) 7.1 9.3 



Lp (g-m"2.yr"M 2.21 0.11 



[P] (mg.m"3) 10 



[chl a] (mg.m"M 1.6 1.4 



S.D. (m) 4.3 5.5 



(yr'M 16 0.24 
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It is not surprising that P concentration should correlate so well 
with the subjective assessment of a lake's trophic status since this 
assessment is normally based on parameters that reflect standing crop 
(e.g. chlorophyll a^ or algal biomass) which are also in units of concen- 
tration. If they had been based on rate processes such as primary production 
then perhaps a P input rate would have been a more appropriate parameter. 

The linkages or sequence of factors affecting P concentration in a 
lake now appear as in Fig. 21 with load being just one of the factors 
controlling P concentrations. This can be expanded into a more detailed 
set of interrelationships (Fig. 22) which is still analagous to Fig. 17 
with its three major groupings. Trophic status, under this new procedure, 
is defined by measuring the parameters of chlorophyll _a, algal biomass, 
Secchi disc and areal hypolimnetic oxygen demand. P concentration is 
represented by the concentration at spring turnover which, in turn, is 
derived from the P loading from a variety of sources and is modified by 
morphometry and hydrology. 

This sequence is the strategy being currently used by the Lakeshore 
Capacity Study within the Ministry of the Environment for arriving at 
answers to the question presented in the introduction. The sequence is 
outlined in the following pages. 

From studies of watersheds in southern Ontario, combined with a 
literature survey, it became apparent that bedrock geology and land use 
were the two major factors controlling P yield or export from watersheds. 
From Dillon's and Kirchner's (1975) results for the two predominant 
bedrock types in Ontario combined with two types of land use (Table 6), 
the disconcerting feature is the broadness of the ranges determined for 
P export. Kirchner (1975) tried to explain some of this variation for 
forested igneous watersheds and found that drainage density (D. = stream 
length/drainage area) explained most of the variation (Fig. 23). 
From this it is possible to make reasonably accurate predictions of the 
P export for a large number of lakes. 

The second source of P is dry or wet precipitation. From Lakeshore 
Capacity Study data and the literature, a range of values of 30 - 70 mg-m ^-yr ^ 
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Fig. 21. Linkages affecting phosphorus concentration 
(m = mass, 1 = length, t = time). 
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Fig. 11. Detailed outline of linkages between watershed and water quality of a lake (from Dillon and Rigler 1975). 
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Table 6: Ranges and mean values for export of total phosphorus from 

43 watersheds. Values include results from present study. 

Results in mg.m"^.yr"^ (from Dillon and Kirchner 1975). 

Land Use Geological Classification 

Igneous Sedimentary 

Forest 

Range 0.7 - 8.8 6.7 - 18.3 

Mean 4.7 11.7 



Forest + Pasture 

Range 5.9 - 16.0 11.1 - 37.0 

Mean 10.2 23.3 



• • 1 ■' 
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Fi<j. 23 Phosphorus export as a function of drainage density (from 
Kirchner 1976). 
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has been established. An averaoe of 50 mc\ P nf-yr"^ is a fairly roliablo 
estimate to use for aeolian loading onto a lake's surface. 

In the past, Dillon and Rigler (1975) estimated the P input from 
human sources by applying 0.8 kg capita"' yr~^ to give a possible annual 
P contribution. This can be qualified, depending on the kind of effluent 
treatment practiced, or in the case of septic tank, tile field systems, 
modified by a soil retention factor (Table 7) as measured by Brandes 
et^ al_. (1974). The problem here is that retention factors for very few 
soil types have been measured. Most apply to artificial mixtures and 
represent only the tile field and filter bed. In fact, determining the 
P input from development in situations where septic tanks and not 
sewage treatment plants exist is one of the purposes of the Lakeshore 
Capacity Study. 

That this problem is complex was demonstrated by Nicholls (1976) 
in a study of Harp and Jerry Lakes on the Precambrian Shield. Harp Lake 
has 73 shoreline cottages while Jerry Lake has none. Both lakes are not 
only morphometrically alike, but also similar in their measurements of 
parameters defining trophic status (chlorophyll a_, Secchi disc, total 
organic nitrogen, areal hypolimnetic oxygen demand, winter oxygen demand) 
(Table 8). The data seem to implicate that the 73 cottages on Harp Lake 
have had little impact on the lake in comparison to Jerry Lake but there 
are qualitative differences that suggest otherwise. Although algal biomasses 
were similar in the two lakes, the blue-green algae were more important 
than diatoms in Harp Lake, while the reverse is true in Jerry Lake. 
Furthermore, the most important diatom in Harp Lake was Tabellaria fenestrate , 
a species that has preceded or accompanied nutrient enrichment problems 
in a number of cases. The same species was entirely absent in Jerry Lake, 
where Rhizoselenia eriensis dominated the community. The real question 
is, "Would both lakes have been the same had there been no cottages on 
Harp Lake?", and that is the question Lakeshore Capacity Study is trying 
to answer in the following manner. 

The Lakeshore Capacity Study has chosen six Precambrian lakes in 
the Muskokas of mid-Ontario with varying amounts of development. Rather 
than express development as the number of cottages, P input to the lakes 
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Table 7: The retention coefficients of total phosphorus from septic 
tile filter beds of different characteristics. Results are 
based on four years of data (from Brandes et al . 1974). 

Filter Bed Rs 



1. 22" Sand 

8" mixture 4% red mud, 96% sand 



7. 15" sand 

15" sand 10^ red mud, 90% sand 

8. 15" sand 

15" mixture, 50% limestone, 50% sand 



10. 15" sand 

15" mixture, 50% clayey silt, 50% sand 



0.76 



2. 30" sand 0.34 

3. 30" sand 0.22 

4. 30" sand 0.48 

5. 30" sand 0.01 

6. 30" sand 0.04 



0.88 



0.73 



9. 39" silty sand 0.63 



0.74 
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Table 8: Comparison of two lakes of similar morphometry, one with cottages 
and one without. (AHOD = areal hypolimnetic oxygen deficit, 
WOD = winter oxygen deficit). 



Harp 



Jerry 



Cottages 



73 



S 



^ax 



^l 



(km') 

(m) 

(lO^mM 

(m) 

(km') 



0.67 


0.50 


12.4 


12.4 


8.3 


6.2 


36 


35 


4.0 


7.3 



[chl a] (mg-m'M 



S.D. 



[P] 



TON 



(m) 



(mg-m"^) 



(mg.l-M 



2.8 
4.4 
13 
0.24 



2.6 
3.9 
16 

0.30 



AHOD 



WOO 



(mg-cm ^-d M 



(mg-cm"'-d"^) 



0.034 
0.022 



0.034 
0.027 



Algae 



Cy.>Ba.^ 



Ba.>Cy.>Cry.* 



* Cy = Cyanophyceae 
Ba = Bacillariophyceae 
Cry = Cryptophyceae 
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was estimated on the basis of the previous information. It has been 
assumed that the per capita P loading for the population is 0.8 kg-yr"^. 
The percent of the total input that could come from people is worked out 
and the progression of selected lakes appears according to Table 9. 
Natural P inputs to these lakes are presently being measured and assuming 
we can reliably predict what a lake's P concentration should be from natural 
loading, then deviations between measured values and the concentrations 
predicted on the basis of the natural inputs must be attributable to 
additional input. Or alternatively, by back-calculating the total loading 
from the lake's total mean P concentration and having measured the natural 
loading, the extra load can be determined by difference and this represents 
the P loading from anthropogenic sources. 

In the preceding discussion it was necessary to convert load into 
a concentration (or vice versa) to provide the necessary translation of a 
lake with trophic status. To accomplish this requires using a P mass 
balance model. The choice of such a model can be made from an increasingly 
complex range of one box models to multi-compartmental models (Figs. 24, 
25, 26). We have used a simple one-box or continuously-stirred tank 
reactor (CSTR) model which only considers input and output rates plus 
sedimentation of a substance, M, in a lake that is instantaneously and 
uniformly mixed (Rainey 1967, O'Connor and Mueller 1970, Vollenweider 
1969, 1975). Other models increase the number of boxes by representing the 
layers of a stratified lake - epilimnion and hypolimnion and even metalimnion 
(Fig. 25) (Sweers 1969). These can be further complicated by adding other 
rate incluencing processes such as production, decomposition, soluble and 
particulate fractions, etc. (Snodgrass and O'Helia 1975, Imboden 1974). 
Then, of course, there are the systems models which try to find solutions 
for coupled sets of differential equations representing the array of 
biological and chemical components, rather than modelling the behaviour of 
a single chemical species (Di Torro et^ al_. 1975, Chen and Orlob 1975). 

The Lakeshore Capacity Study model can be described by the equation 
in Fig. 27 where M^ is the outflow concentration (Dillon and Rigler 1975). 
Mq was originally equated with the average lake concentration but it can be 
approximated by the second expression (Chapra and Tarapchuk 1976). This 
empirically states that mean annual outflow total phosphorus, [P] is 905K 
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Table 9: Estimated portion of total phosphorus supplied (J,) by 

anthropogenic sources (J.) to Lakeshore Capacity Study Lakes 
(0 = P supply in mg yr~'). 



Lake 'ip/Jj{%) 



Jerry 



Red Chalk 



Chub 26 



Blue Chalk 27 



Dickie 54 



Harp 61 



CLASS 2 - STRATIFICATION 



2q) Conservative 
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WINTER 





^'TM 



eg. Sweers - pollutants in Greot Lakes 
2b) Non- conservative 




^M _J~M 



o[m] 



'TCm] 



Fig. 25 Schematic representation of mass balance models that consider 
thermal stratification. 
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CLASS I ONE -BOX MODELS 



la) Conservative 



L. 

2 




^<rM 



eg. Roiney - pollutants in Great Lakes 
O'Connor a Mueller - CI in Great Lakes 



lb) Non- conservative 

_L 

z 




^^M] 



eg. Vollenweider - P in lakes 



Fig. 24 Schematic representation of "one-box" mass balance models. 



-44- 



CLASS 3 - IN- LAKE EVENTS INCLUDED 



[M]^M 



SETTLING 



DIFFUSION 



\r 






settling 
_jL__ 



-4 



[M]?±[Mg 



M? 



eg. Snodgrass \ „ :„ i«i,«e 

CLASS 4 - SPATIAL INHOMOGENEITIES S 

BIOLOGICAL COMPARTMENTS 



eg. Thomann, DiToro 



Fig. 26 Schematic representation of mass balance models that include 
in-lake mechanisms. 
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d[M] = J^ - cr[M] -_Q_ [m]o 
dt 2 V 

[P]o = 0-9 [P] 

[P]s8= l/z(<r+0-9^) 

or 

[P]^^= L(|-R)/o-92^ 



Fig. 27 Mass balance equation for a substance M, where L is the loading 
rate, z the mean depth, o the sedimentation rate, Q the outflow 
volume, V the lake volume, p = Q/V, [M]© the outflow concentra- 
tion, R the retention coefficient, [PJss ^^^ steady-state 
phosphorus concentration (where M e P). 
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of the value of the mean lake total [P]. The steady-state can take on two 
forms, the second one containing the parameter R or the retention 
coefficient (R = the fraction of the total P input that does not leave a 
lake = a/(a+p)) (Dillon and Rigler 1974a). 

The model was first tested by measuring all the required parameters 
(L, z, p, R) for a series of lakes excluding the 0.9 factor (Dillon and 
Rigler 1974a). The actual measurements compared very well with the 
predicted P concentration at spring turnover (a mean deviation of 1.5 yg 1"^) 
(Table 10), but when the 0.9 factor is included the predicted results are 
improved (a mean deviation of 1.0 yg I'M- Literature cases have also 
yielded excellent results using the steady-state equation. 

For this or any other mass balance equation to be of general use, 
it should be possible to predict as many of the input variables as possible. 
Methodology for estimating P loading (L) was previously outlined and z 
is usually obtainable for any lake using morphometric maps. This still 
leaves two parameters to estimate, p and R. Lake flushing rate or p can 
be approximated by deriving a water budget for a lake from precipitation, 
evaporation and evapotranspi ration data. (Dillon and Rigler 1975). Such 
data have usually been isoplethed using long term averages for various 
geographic regions by government agencies and are readily available. 

R has also been related to hydrologic parameters but with moderate 
success such as in Fig. 28 of R versus p times i (p'Z = hydraulic loading = 
q's) (Kirchner and Dillon 1975). Larsen and Mercier (1976) have tried 
other combinations of these parameters but none of these appears to have 
strengthened predictions of R. Before more empirical modelling is attempted 
in this area, some fundamental questions have to be answered, such as 
"Does R change as L changes?". 

The preceding paragraphs have outlined a basic technique that 
allows prediction of the mean vernal P concentration in lakes. What 
remains to be done is to link this P concentration to factors indicative 
of trophic status. 

One of the first relationships, studied by Dillon and Rigler (1974b) 
following Sakamoto's lead (1966), was the relationship between P concentration 
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Table 10: Comparison of measured phosphorus concentration at spring 
overturn with that predicted by a simple mass balance model 
(see Dillon and Rigler 1974a). 



[P] (mg.m"^) L(l-R)/2p 



Cameron 12.4 10.8 



Four Mile 12.0 8.3 



Bob 8.S 7.1 



Boshkung 6.1 5.4 



Halls 4.3 4.0 



Beech 8.5 7.8 



Maple 9.3 7.4 



Pine 10.5 8.4 



Cranberry 9.9 6.9 



Moose 7.1 6.8 



Hall burton S.3 6.1 



|a| 1.5 




, v=IOmyr-i • 



J u 



10 20 40 60 



80 100 120 140 160 180 200 220 240 

Qs 



Fig. 28. Relationship between phosphorus retention coefficient (Rp) and hydraulic load (qs = zp}. The settling 
velocity v can be shown to be equal to zp. 
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at spring turnover and mean sumer chlorophyll £. Fig. 29 is the original 
data set used, including Sakamoto's data (Dillon and Rigler 1974b), while 
Fig. 30 shows the addition of later points. Although the relationship 
appears to work very well, it is a log-log relationship, so it has limited 
prediction capabilities at both of its tail-ends. 

Because of the relative ease and effort required to measure 
chlorophyll £, many P-chlorophyll applets have been published in the 
literature (Fig. 31) and great differences between one another have been 
shown. Nicholls and Dillon (1978) have attempted to analyze the reasons for 
such differences (see Appendix A for further discussion). 

Secchi disc depths have been estimated from chlorophyll £ (Fig. 32) 
by a number of people. There are three factors affecting light extinction 
and thus the correlation between Secchi disc and chlorophyll a: 1) particu- 
late material containing plant pigments, 2) particulate material not containing 
plant pigments and 3) dissolved colour. The second factor may be of little 
importance in Precambrian lakes, but colour frequently is. Colour can be 
factored out quite easily by multiple regression analysis. 

Relationships between average total phosphorus concentration and phyto- 
plankton biomass (cell volume) and between phytoplankton biomass and Secchi 
disc visibility are generally better defined than those substituting chlorophyll 
as an indicator of phytoplankton biomass. However, chlorophyll will surely 
remain a popular measure of algal biomass due to the economics and time 
constraints of most research programnes. 

Ideally, one would like to measure parameters such as P loading and 
concentration and predict the genera or species of zooplankton, phytoplankton 
and macrophytes (etc.) to be found in a lake. Nicholls (1977) found that 
the greater the algal biomass in a lake, the less important the Chrysophyceae 
are as contributors to this biomass (Fig. 33). When this is coupled with 
P data, it is apparent that the Chrysophyceae disappear from a lake ecosystem 
at a P concentration greater than 30 ug 1"^ (Fig. 34) for reasons 
explained earlier. This kind of information is of particular value to 
water managers. For example, the small motile Chrysophyceae are an important 
nutritive source for the zooplankton. In fact, their small size and hence 
rapid turnover suggests their importance in the production ecology of a lake 
is probably much greater than measurements of their biomass would indicate. 
The loss of this group from a eutrophying lake clearly has important 
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Relationship between total phosphorus concentration at spring 
overturn and average summer chlorophyll a concentration {from 
Dillon and Rigler 1974b). 
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30 Additional data shown with the original line generated in Fig. 29 
(Dillon and Rigler 1975b). 
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Fig. 31. Phosphorus - chlorophyll curves reported in the litere-.jre (Nicholls and Dillon 1978) 
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Relationship between Secchi disc depth and chlorophyll a^ 
concentration (summer average) in the euphotic zone of 
Ontario lakes. Data are taken from the Ontario Ministry 
of the Environment recreational lakes progranme. 
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Ftg 33 Relationship between total ph^'toplankton biomass and percentage of the biomass contributed by Chrysophyceae 
(from Nicholls 1977). 



CD 



30- 



o 

CD 

g20- 

o 



a 
o 

CO 



10- 



o 



0- 




Fig. 34. Relationship between total phosphorus content of lakes 
and the proportion of the algal biomass contributed 
by the Chrysophyceae. 
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implications on the food chain and ultimately, the fisheries. Consequently, 
Nicholls (1977) has developed a nomogram (which also considers the relative 
importance of the Chrysophyceae, expressed as a percentage of the total 
phosphorus) to interrelate these indices of trophic status for predictive 
purposes (Fig. 35). An example of the use of these techniques is 
outlined. The P budget for Gravenhurst Bay was estimated before and after 
P removal (Table 11). The anthropogenic load from the sewage treatment 
plant was measured while the other figures were estimated or extrapolated. 
The 1975 predicted and measured values for several "trophic status indices" 
(Table 12) suggest that changes in P concentration can yield quantitatively 
predictable results. 

The inherent assumption in the use of a phosphorus-phytoplankton 
relationship for predictive purposes is that the total phosphorus measure- 
ment includes all phosphorus in solution and suspension which, if not 
already incorporated into algal cells, is potentially available. Clearly 
there may be some forms of phosphorus in lake water which are not available, 
yet it does not seem practical to determine routinely the "biological P" 
in lake water since such studies involve rather sophisticated bioassay 
techniques. It is nevertheless important to recognize in nutrient budget 
studies that a portion of the phosphorus load to a lake may never be available 
for plant growth in the lake (for example see our results of algal assays 
on Lake Timiskaming inflowing streams) (Fig. 36). 

Although there is now much data to provide strong evidence that 
algal growth is controlled by phosphorus concentration in most lakes, 
especially those of the Precambrian Shield, there is also evidence that 
very productive waters may be nitrogen controlled (e.g. Nicholls 1976a). 
Additionally, we have completed laboratory algal assays on the inner Bay 
of Quinte and Thompson Lake and have demonstrated nitrogen limitation in 
these highly eutrophic waters (Fig. 37, 38 and 39). 

These findings suggest important implications for the recovery of the 
Bay of Quinte following implementation of controls on phosphorus loading 
from municipalities within the Quinte basin. Because the inner bay 
photoplankton is nitrogen limited and the outer bay phytoplankton, phosphorus 
limited, the response to decreased phosphorus loading from sewage treatment 
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Fig. 35 Nomogram facilitating intercomparisons of the most reliable trophic status indicators. To illustrate the 
use of the nomogram, for predictive purposes, Balsam Lake average total phosphorus concentration and the 
resulting predicted values of total phytoplankton biomass, Chrysophj'ceae abundance and Secchi disc 
visibility are indicated. The arrows indicate the measured values. 
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Table 11: Estimated Phosphorus budget for Gravenhurst Bay (from 
Dillon et al . 1978). 



Source 



1971 



1975 



Land drainage (kg) 



180 



180 



Cottages (kg) 



520-1290 



280-690 



Municipal (kg) 



2580 



680 



Precipitation (kg) 



150-310 



150-310 



Total (kg) 



3430-4360 



1280-1860 



Load (g-m"^.yr ^) 



0.82-1.06 



0.31-0.45 
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Table 12: Comparison of predicted and measured trophic status parameters 
for Gravenhurst Bay following P removal. 



Measured Predicted 



[P]sp (mg-m"3) 26 20-29 



[PJ55 {mg-m-3) 20 17-25 



[chl a] (mg-m"^) 5.0 5.6 



S.D. (m) 3.9 2.7 (3.0 from nomogram) 



AHOD (mg-cm-2.d"M 0.025 0.035 



Phytoplankton /^3.1-M 05 2 

biomass inm ' J ^.t) ^.u 



Chrysophyceae {% of total) 11 15 
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Fig. 36. 



Laboratory algal assay of stream water inputs to Lake Timiskaming. . 
Keyed as follows: A, Alligator Creek; W, Waki Creek; L.M., Lower 
Moffat Creek; L.O., Lower Otterskin Creek; B, Blanch River; 
Q, Quince River; U.M., Upper Moffat Creek; U.O., Upper Otterskin 
Creek; T. Lake Timiskaming (mid-lake). To some extent, the poor 
fit of the data suggest that a portion of the total phosphorus 
measured at some stream locations is not available for algal 
growth. 
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Growth of Monoraphidium griffithii in filtered water from the 
inner Bay of Quinte (control) and showing the increased yield 
of the water with added NaNOs (1,125 mg N/1). Three additional 
samples spiked with KH2PO., (50 pg P/1 ) showed no significantly 
increased yield over the control. 
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Fig. 38. Growth of Monoraphidium griffithli in filtered water from the 
outer Bay of Quinte (control) and shovnng the increased yield 
of the water with added KH2PO1, (50 ug p/1 ) . Three additional 
samples spiked with NaNOa (1,125 mg N/1) yielded lower maximum 
biomass than either the control or the P-spiked samples. 
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Growth of Selenastrum caprice rnutum in filtered water from Thompson 
Lake showing growth response to additions of inorganic phosphorus 
and nitrogen. A similar experiment showed no response to additions 
of iron and manganese. 
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plants should be noticed first at the outer bay locations. In the inner 
Bay of Quinte, phosphorus will have to be established as the controlling 
nutrient before further reductions in concentration permit declines in 
algal biomass. 

Nitrogen limitation of phytoplankton in highly eutrophic lakes 
generally does not curtail the obnoxious bloom- forming species 
(e.g. Aphanizomenon and Anabaena spp.) which can fix gaseous nitrogen. 
In fact, in Thompson Lake (15 km N of Toronto), these bloom- forming species 
were dominant during most of the summer period prior to attempts to 
control them by artificially destratifying the lake (Fig. 40). By means 
of their gas vesicles, blue-green algae are unique in their ability to 
alter their buoyancy and select depths which are optimal for their growth, 
A stable water column is therefore essential for the effective use of their 
buoyancy regulating mechanism. That certain blue-green algae can be 
essentially eliminated from a thermally stratified lake by mixing the water 
has been known for some time (e.g. Sirenko et^ al_. 1972, Karasik and 
Stovbun 1974). 

When a eutrophic lake is artificially mixed (for example, by 
aeration), blue-green algal cells are forced to spend a disproportionate 
amount of time in the deeper, "aphotic" water which promotes excessive 
formation of gas vesicles within the cell. The buoyant cells are then more 
prone to surface floatation where they are killed by sunlight (ultraviolet 
exposure and/or photo-oxidation of pigments). During 1974, when Thompson 
Lake (Aq = 3.9 ha, Z = 26 m) was being adequately mixed by the aerator, 
blue-green algae were almost non-existent in the lake (Fig. 41). Instead 
the lake was populated by small coccoid green algae, cryptomonads and 
diatoms, forms which fit well into aquatic food chains and promote 
development of larger fish-food organisms. Owing mainly to their large 
colonial habit and the presence of a protective mucillagenous sheath, most 
blue-green algae are not an important component of the food chain. 

It is worth emphasizing that lake systems are rarely so simplified 

that a single stress illicits a single response, and in some eutrophic 

lakes, lake response to artificial mixing apparently can be much more 

complex. Heart Lake (A„ = 14.3 ha, Z^^^ = 10.4 m), which has also had a 
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Fig. 40. Phytoplankton of Thompson Lake during 1972. Note the high biomass of blue-green algae (Cyanophyceae) 
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F1g. 41. Phytoplankton of Thompson Lake during 1974, Note the virtual absence of blue-green algae (Cyanophyceae) 
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history of blue-green algal problems, has been aerated sinco 1975. 
Artificial destratification again apparently resulted initially in replacement 
of blue-green algae by "food-chain functional" forms such as small greens 
and cryptomonads. However, an additional effect of aeration was to increase 
the habitat of Daphnia pulex by replenishing the oxygen content of the 
deeper waters (Strus 1976). Because the D^. pulex could now take refuge 
in the dark, deeper, oxygenated water from planktivorous fish, they were 
therefore able to thrive under the conditions created by aeration. The 
development of the Daphnia population was so successful that the desirable 
forms of algae were consumed at a rate almost equal to their production. 
By this stage in the sequence of events, the large dinoflagellate 
Ceratium hirundinella had appeared {as it usually does in July and August 
in eutrophic lakes). It is too large to be grazed by Daphnia pulex , 
and in the absence of other algae, C^. hirundinella experienced no competition 
for light and nutrients and a massive population developed (Fig. 42). The 
sudden death of this massive dinoflagellate population resulted in a 
sumrnerkill of the lake. 

Clearly, some lake mangement techniques have complex ramifications 
and a thorough understanding of lake processes in a lake-by-lake assessment 
would seem to be a necessary requisite to any successful management 
programne. 

In summary, the points that we have tried to make are the 
following: 

1) That P is the element controlling trophic status in most lakes and 
that biological changes can be induced by changing P concentration by 
even very small amounts. Under some highly artificial conditions, N or 
C may control algal growth. N limitation especially can be easily 
detected with a simple laboratory algal assay. 

2) That the input of P to a lake can be estimated without a lot of 
collection of raw data. 

3) That lake morphometry, hydrology and phosphorus load can be combined 
in a mass balance model to give a good prediction of phosphorus 
concentration. 
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Fig. 42. Seasonal distribution of total phytoplankton and of Ceratium hirundinella in Heart Lake during the first 
and second years of artificial destratification. 
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4) That P can be used to let us estimate pliytoplankton bioiiiass and some 
characteristics of the composition which liave important implications 
for lake management (e.g. hypolimnetic dissolved oxygen depletion 
rates, efficiency within the food v/eb of the lake, water clarity, etc.) 
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APPENDIX A 

A literature search has yielded 16 different P-chlorophyll 
relationships. The following is a condensed version of Nicholl's and 
Dillon's (1978) explanation of this variation, for which they concluded 
that there are two basic reasons: 

1 . Differences in laboratory and field techniques 

The upper line in Fig. 31 is Sakamoto's original plot (1966) and his 
analyses are of total chlorophyll rather than just chlorophyll a^. The 
four plots grouped in the middle used spring P levels while almost 
all others used mean summer or annual P concentrations. This difference 
clearly explains some of the variation in the relationships, particularly 
in view of the fact that Chapra and Tarapchak (1976) found that summer 
total phosphorus concentrations were consistently lowe'^ than spring 
values by a factor of 0.9 (Fig. 43), an observation that is well supported 
by Ontario lakes (Fig. 44). 

In short, no two of the 16 studies used the same sampling periods, 
laboratory procedures and the same sampling zones in the lakes studied. 

2. Chlorophyll a^ is no more than a crude approximation of algal biomass 

It is now known that chlorophyll a_ content per unit of algal cell volume 
varies from 0.1 to 9.7%. Fig. 45 is a plot of results from a study 
of the Kawartha Lakes in southern Ontario (Nicholls 1977). It is 
obvious that switching from chlorophyll a_ concentration to volumetric 
algal biomass provides a stronger relationship. Similarly, Secchi disc 
biomass plots are better correlated than Secchi disc - chlorophyll a_ 
plots (Fig. 46). This superiority of 'measured' biomass is to be expected 
and is the better measurement of biomass even though a simple chemical 
measurement, like chlorophyll a_ or ATP, will still have a place in 
limnological studies because of their ease in execution. 
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Fig. 43. Relationship between total phosphorus concentration measured 
at spring overturn and that calculated on an annual basis 
(from Chapra and Tarapchak 1976). 
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Relationship between total phosphorus concentration measured at spring overturn and that calculated on 
an annual basis for lakes in southern Ontario. 
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Fig. 45. Chlorophyll a concentration and phytoplankton biomass vs total phosphorus content for the Kawartha Lakes 
{Nicholls 1977). 
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Fig. 46. Secchi disc depth vs chlorophyll a concentration and phytoplankton biomass for the Kawartha Lakes 
(Nicholls 1977). 
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